Interaction-driven topological insulator states in strained graphene 
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The electronic properties of graphene can be manipulated via mechanical deformations, which 
opens prospects for both studying the Dirac fermions in new regimes and for new device applications. 
Certain natural configurations of strain generate large nearly uniform pseudo-magnetic fields, which 
have opposite signs in the two valleys, and give rise to flat spin- and valley-degenerate pseudo Landau 
levels (PLLs). Here we consider the effect of the Coulomb interactions in strained graphene with 
uniform pseudo-magnetic field. We show that the spin/valley degeneracies of the PLLs get lifted 
by the interactions, giving rise to topological insulator-like states. In particular, when a nonzero 
PLL is quarter- or three-quarter filled, an anomalous quantum Hall state spontaneously breaking 
time-reversal symmetry emerges. At half-filled PLL, weak spin-orbital interaction stabilizes time- 
reversal-symmetric quantum spin-Hall state. These many-body states are characterized by the 
quantized conductance and persist to a high temperature scale set by the Coulomb interactions, 
which we estimate to be a few hundreds Kelvin at moderate strain values. At fractional fillings, 
fractional quantum Hall states breaking valley symmetry emerge. These results suggest a new route 
to realizing robust topological insulator states in mesoscopic graphene. 

PACS numbers: 73.43.-f, 72.80.Vp 



Introduction. Graphene provides a unique exam- 
ple of a 2D membrane hosting a high-mobility two- 
dimensional electron gas Mechanical deformations 
of the membrane strongly couple to the electron motion, 
generating local potentials and gauge fields . The 
interplay of the mechanical and electronic degrees of free- 
dom have been mostly discussed in relation to the trans- 
port properties of graphene samples which are often lo- 
cally strained and rippled T|. More recently, it was real- 
ized that strain engineering provides a promising way to 
control the band structure and the spectrum of graphene, 
as well as to create new devices [7H9|. 

Non-uniform strains locally shift the positions of the 
K, K' Dirac nodes in the opposite directions; this is 
equivalent to the local vector potential, and, in general 
creates local pseudo-magnetic field, which has opposite 
signs in the two valleys This fact implies that the 

time-reversal symmetry (TRS) is preserved by strain, as 
well as that translational symmetry is broken, since the 
states in two valleys generally have different magnetic 
translation operators. Recently, it was shown that cer- 
tain natural strain configurations lead to a large, nearly 
uniform pseudo-magnetic field [13, HI- This drastically 
modifies the spectrum of graphene, creating a sequence of 
quantized, fourfold (valley and spin) degenerate pseudo 
Landau levels (PLLs). Recently, PLLs were observed in 
locally strained graphene [13, [13] ; the effective pseudo- 
magnetic field was as high as 300 T [l2j . 

Main results. We study the effects of the elec- 
tron interactions on the physics of partially filled PLLs 
in strained graphene. We show that the approximate 
four-fold valley and spin degeneracy of the PLLs gets 
lifted, giving rise to incompressible topologically non- 
trivial states. We find that in n 7^ PLLs, where n 



numbers the PLLs in the standard way m , three distinct 
types of states emerge - anomalous quantum Hall (AQH) 
state [l3 | breaking TRS, TRS-preserving quantum spin- 
Hall state (QSH) [l5j, as well as a trivial insulating state. 

The AQH and QSH states are characterized by chiral 
and counter-propagating edge states, respectively. The 
edge states are insensitive to the crystallographic edge 
type; this should be contrasted with the case of incom- 
pressible states corresponding to the complete filling of a 
PLL, which, in general, have gapless edge states only for 
particular edge orientations [10|. 

The mechanism leading to the formation of the broken- 
symmetry states is reminiscent of the quantum Hall ferro- 
magnetism (QHFM), for the case of graphene discussed 
in Ref. 16[. Similarly to QHFM, in our problem val- 



ley/spin split states are favored by the exchange interac- 
tions; however, owing to the different underlying symme- 
tries, the properties of the resulting split states are very 
different from the usual QHFM states. In the leading 
order the Coulomb Hamiltonian of a partially filled PLL 
has a Z2 X SU{2) symmetry; Z2 is the time reversal, which 
interchanges valley and spin species, while SU (2) corre- 
sponds to the spin rotation symmetry. This should be 
contrasted with the case of an SU (4)-symmetric QHFM, 
where any linear combination of the two valleys can be 
occupied [l6|). 

Below we focus on the case of non-zero PLL (the case 
n = 0, in general, does not give rise to topologically non- 
trivial states, as discussed at the end of the paper; also, 
see Ref. At partial filling / = 1, 3 of the n 7^ 

PLL quartet, the ground state is an Ising-type valley fer- 
romagnet. The strong Ising-type anisotropy stems from 
the lack of the exchange interaction between the states 
in K and K' valleys, which results from the opposite chi- 
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rality of Landau orbitals in K and K' valleys. The valley 
ferromagnet at / = 1,3 breaks Z2 TRS, and is charac- 
terized by the quantized Hall conductivity a^y — ±e^//i, 
with the +(— ) sign for the K{K') valley polarization. 

At partial fiUing / = 2, a QSH state is stabilized by 
weak SO interactions. The main effect of the SO in- 
teraction is to make the wave function orbitals of one 
pair of PLL related by TRS different from those of the 
other pair. Because of that, one of the pairs of PLLs 
has a lower exchange energy; the ground state at / = 2 
thus corresponds to filling that pair of PLL. Such a state 
can be qualitatively viewed as two filled LLs for spin- 
up and spin-down electrons, which experience opposite 
magnetic field; this state preserves TRS and is charac- 
terized by the quantized spin-Hall conductance and zero 
charge-Hall conductance. In order for the QSH state to 
be energetically favorable, the SO-induced exchange en- 
ergy anisotropy should exceed weak anisotropy due to the 
lattice effects, which favors ferromagnetic spin-polarized 
state, and which is estimated below. 

The distinctive property of the interaction-induced 
QSH and AQH states is their robustness. They per- 
sist up to a temperature scale set by the Coulomb in- 
teraction, which, as we estimate below, can reach 200K 
at a moderate strain-induced pseudomagnetic field of 
B = lOT This is much larger than the typical tem- 
perature scales set by bare SO interaction, which char- 
acterize the QSH states in graphene 15| and in HgTe- 
based quantum wells [18(, JJiJ ; for example, the latter per- 
sists only to temperatures of several Kelvin. We expect 
that the robustness of the interaction-driven QSH states 
will be advantageous for their manipulation and for the 
studies of their edge states properties. Potentially, QSH 
states in strained graphene could also be used in spin- 
tronics applications. 

Model. The general single-particle Hamiltonian of the 
strained graphene is given by: 

Hk{K') = voi'HPxO-xTz + Pya-y)i>, Pi ^ Pi ± sA^, (1) 

where vq is the Fermi velocity, and K{K') correspond 
to the two valleys of the Dirac fermions, located in the 
opposite corners of the Brillouin zone; cr^, are the Pauli 
matrices acting in the sublattice and valley spaces. The 
vector-potential is related to the strain gradient • The 
above Hamiltonian does not contain spin-orbital terms, 
which will be included below. 

We will be interested in the strain configurations that 
give rise to nearly uniform pseudo-magnetic field; such 
configurations were described in detail in Refs. [lol. 
The corresponding vector potential is given by 



by fl: 



= -^By, Ay = ^Bx. 



(2) 



Such pseudo-magnetic field breaks the spectrum of 
graphene into a sequence of PLLs with energies given 



En = sgn{n)y/\n\eo, Eq 



V2h 



(3) 



where ^ = y ^ is the effective magnetic length set by 
B, and vq is is assumed to be constant, as appropriate 
for strains which are not too strong. In what follows, 
without loss of generality, we take B > for definiteness. 

In our analysis of the Coulomb interactions, we will 
need PLL wave functions. The details can be found in 
the Appendix; here we will just use the final result. The 
n 7^ PLL wave functions of the Hamiltonian in K valley 
are given by: 

■ipn,ni(z) = (</'|«|-l,m(2), Sgn(n).^|„| ^„ (z)) , (4) 

where ipn,m are the wave functions in the nth non- 
relativistic LL (see Appendix for the definition). The 
zeroth LL wave functions have a different form, with only 
bottom spinor component being non-zero: 



^Pn,miz) = (0, (po,miz)) ■ 



(5) 



The wave functions in the K' valley have the same form 
as those in the K valley ((24)) (notice that this is different 
from the case of graphene in real magnetic field, where 
one would have to interchange the upper and lower spinor 
components). However, notice that the magnetic oscilla- 
tor wave functions ipn,m describe the opposite direction 
of the cyclotron motion compared to the K valley. 

The Coulomb Hamiltonian. In our analysis of the 
interaction effects, we will neglect the effects of Landau 
levels mixing. This amounts to projecting the density 
operators onto the partially filled rtth LL. The effective 
Hamiltonian of the partially filled PLL can be written 
down in analogy with the case of real magnetic field (for 
a review, see Ref. [20l|): 

H ^ ^Y.^{<l)p{ci)p{-ci), p(q) =5]pK.(q), (6) 

KS 

where k = K, K' , s =t, i label valley and spin states, and 
PKs(q) are the density operators projected onto the nth 
LL. They have a different form in the K and K' valleys: 



F„(q)G, 



(7) 



m.ni' 

(8) 

where we defined q,q = Qx HQy, being measured in units 
of and G is given by 

G„,„(9) = (m!/rn'!)i/2(_,^/V2)™'-™L™'-™(|(7|V2), 
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being the generalized Laguerre polynomial; i^„(q) are 
the graphene form-factors, which encode the structure 
of PLL wave functions, and are identical to the case of 
magnetic field [3| 

K(q) = i [L|„|_i((?V2) +i|„|(gV2)] . (9) 

The different structure of the projected density operators 
(|7l8p reflects the fact that the pseudo-magnetic field has 
the opposite direction in the two valleys. This difference, 
although it may appear minor, has important implica- 
tions for the energetics of the broken-symmetry states. 

Competing phases. Similarly to the case of QHFM, 
the physics of the partially filled PLL is governed by the 
exchange interactions. Below we show that the exchange 
energy is minimized when / sub- levels in K and/or K' 
valley and with arbitrary spin projection are filled (that 
is, we are dealing with an Ising-type valley ferromag- 
net, and no linear combinations of the two valleys are 
allowed) . This stems from the different cyclotron motion 
direction of the wave functions in K, K', which causes 
the exchange contribution due to the correlations be- 
tween different valleys to vanish [21]. We also show that 
the ground state degeneracy is lifted by weak SO inter- 
actions, which break SU{2) symmetry of the Coulomb 
Hamiltonian. This, generally, chooses a particular spin 
orientation at / = 1,3, and favors QSH state at / = 2. 

Partial filling f = 1. We will start our analysis with 
the case of partial filling / = 1. As we shall see, some 
of the results obtained here will be also useful for the 
case f — 2. By analogy with the case of the quantum 
Hall ferromagnet, let us consider the following trial wave 
function: 

i*>=n4jf^), (10) 

where \Q) is the vacuum, and d is obtained by a unitary 
rotation U of different valley and spin states: 

= y ] ^lisCj^s.nm- (H) 

KS 

Evaluating the average of the Hamiltonian (jG)) over the 
state (fTO|) . and noting that the exchange integrals van- 
ish between states in different valleys in a macroscopic 
sample, we find the exchange energy per particle in the 
relevant PLL [2l[: 

E ^ -A„ [{riKt + TiKif + {riK'^ + riK'if] , (12) 

where n^s = UksU^s is the filling factor of the ks PLL, 
and 

From the Eg. (fT2l) . we see that the lowest energy states 
are completely valley polarized, uk^ + uki = 1 , riK^i + 
n-K'i = (or vice versa). 



The ground state is a valley ferromagnet that breaks 
TRS. Such a state is characterized by the quantized Hall 
conductivity and chiral edge states. The domain walls 
between domains with opposite valley polarization will 
also carry two chiral ID modes. 

What lifts the spin degeneracy of the ground state? 
The spin-orbital interactions, which entangle orbital and 
spin degrees of freedom, do not always lift the single- 
particle degeneracy of the Landau levels. However, they 
do modify the PLL wave functions for different spin pro- 
jections (TRS is still preserved, such that there are two 
pairs of identical orbitals) . This modifies the correspond- 
ing form-factors, breaking the SU(2) invariance of the 
energy functional, and favoring some particular spin pro- 
jection. 

We illustrate this general effect by considering the in- 
trinsic SO interaction [l5| : 

Hso = AsoTzCTzSz, (13) 

where Si acts in the spin space. This type of SO inter- 
action conserves Sz component of spin, and therefore is 
the easiest to analyze; other types of SO have the same 
qualitative effect. 

The intrinsic SO (|13p modifies the wave functions in 
the K valley ([M)) . making the spin- up and spin-down 
wave functions different: 

V'ra,m(^) ~ (cos (z) sin ^„ (z 

)) , (14) 

V'im(^) = (sin^V|«|-l,m(^) COS^.^|„|^„(z)) , (15) 

where tan^ = sgn{n)y^l + (Aso/eo)^ - ^so/so- Simi- 
larly, the wave functions in the K' valley will be modified, 
in such a way that the TRS is preserved: K' f {K' \.) 
state wave function will have the same form as the 
K I {K t) wave function. Thus, there are now two pairs 
with different spatial wave functions. 

The difference in the wave functions makes the form- 
factors of the two spin species different, leading to the 
difference in the exchange energy. The new form-factors 
are given by: 

F,f^(q) = F5^(q) = cos^ eL|„|_i(gV2) W ei|„| (gV2), 

(16) 

F,'^H^) = FfHd) ^ sin2eL|„|_i(gV2)Wei|„|(gV2). 

(17) 

This breaks the spin-rotational symmetry of the 
Coulomb Hamiltonian, generating an analogue of the 
Zeeman interaction in the exchange energy functional: 

E = -A„(nK^ + riKif - S^{nK^ - uki), (18) 

where we assumed that the electrons have already been 
valley-polarized in K direction by their strong Coulomb 
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repulsion, and the Zeeman-like term is given by, to the 
leading order in Aso/eq, 



A 



so 



^y(g)e-V (if„|-i(«V2) - 4.|(<zV2)) 

(19) 

For \n\ = 1, evaluating the integral, we obtain Sz = 
, where a = e^/hvoe is the coupling constant. 
This anisotropy, although small, fixes the spin polariza- 
tion to be = 1. Alternatively, a state with the same 
energy could be obtained by filling K' I PLL. 

Partial filling f — 2. Now we proceed to the case 
/ = 2. In the leading order (neglecting SO interaction), 
there is a manifold of degenerate ground states: any two 
non- equivalent PLLs {ks, k's', where k, k' — K,K', and 
s, s' are linear combinations of spin- up and spin-down 
species) can be filled. Similarly to the case f = 1, this 
degeneracy is lifted by the SO interactions, which modify 
the PLL wave functions and reduce the symmetry of the 
Coulomb Hamiltonian. For the case of intrinsic SO, the 
exchange energy is minimized when a time-reversed pair 
of PLL, K t, K' i, is filled. 

Such a state is TRS-preserving, and can be viewed as 
a combination of filled spin-up and spin-down LLs sub- 
ject to an opposite magnetic field. Via a Laughlin-type 
argument, we conclude that such a state is characterized 
by two counter-propagating spin- filtered edge states , 
and a quantized spin-Hall conductivity, asH — o-} 



xy 

2e^lh. Moreover, the gapless edge states are 
protected as long as TRS is preserved; this can be seen 
by employing Kane and Mele's S'-matrix argument [l5| . 
This argument shows that the TRS forces the off-diagonal 
elements of the edge states S'-matrix to be zero, and thus, 
backscattering is prohibited as long as TRS is not broken. 

Case of Rashba interaction. Other types of SO inter- 
action (for concreteness, we focus on the Rashba term), 
do not conserve Sz projection. Despite this, their effect 
on the / = 2 is similar to the s^-conserving intrinsic 
SO: they break the spin rotational symmetry, making 
the wave functions of two time-reversed PLL pairs dif- 
ferent [111 . One of the pairs of PLLs has lower exchange 
energy, and the ground state corresponds to filling that 
pair. 

Importantly, such a state is also a quantum spin-Hall 
insulator, with a pair of protected gapless edge states. 
This is because it can be adiabatically connected to the 
quantum spin-Hall insulator in the presence of only in- 
trinsic SO interaction (when SO interaction is adiabat- 
ically changed from pure intrinsic type to the Rashba 
type, the exchange gap does not close). This should 
be contrasted with the Kane-Mele model [11] , where the 
Rashba interaction weakens and eventually destroys the 
quantum spin-Hall state. 

Tc estimate. The critical temperature at which the 
interaction-induced topological states set it, is deter- 
mined by the exchange gap, Tc ~ A„. For n = 1 PLL 



we obtain A„ = 55 v f fl' '^^^re we restored the dimen- 
sional factors. Taking e « 5 (the intrinsic screening of 
graphene [J) gives A„ « 200K at S = 10 T. This scale 
far exceeds the bare SO interactions Ago; thus, we ex- 
pect the interaction-induced QSH state discussed above 
to be significantly more robust than the QSH insulator 
in non- interacting systems 15 , . Such a state will also 
show greater robustness with respect to disorder. 

Concluding remarks. Before we conclude, several 
comments are in order. First, much of the analysis above 
carries over to the splitting of the n = PLL, previously 



considered in Ref. [17|. At partial filling / = 1,3 of the 
zeroth PLL, we expect valley-polarized states, while at 
/ = 2 weak SO interactions would favor TR-symmetric 
state (also, note that intrinsic SO splits n = PLL al- 
ready on the single-particle level). However, in contrast 
to the case of n ^ PLLs, we do not expect these states 
to exhibit protected edge states; this is because, in gen- 
eral, zeroth PLLs in K{K') valleys cannot be character- 
ized by a quantized Hall conductivity. 

We also expect the electron interactions to give rise to 
fractional quantum Hall states, which break TRS sponta- 
neously. In particular, we expect robust fractional states 
to emerge at partial fillings / — fc + 1/3, / — k + 2/3 , 
k = 0,1,2,3, of the n = and n = 1 PLLs [IH, 
while higher PLLs form-factors favor charge-density- 
wave states. Although a possibility of TRS-preserving 
fractional states [26i | is intriguing, within the variational 
approach for the case of pure long-range Coulomb inter- 
action we found no evidence that such states could be 
energetically favorable 2l|. Notice that TRS-preserving 



fractional states found in [17| required carefully chosen 



combination of lattice-scale interaction parameters. 

Discussion. In summary, we have shown that nonzero 
pseudo Landau levels in strained graphene host a number 
of broken-symmetry states, including quantum Hall-like 
states which occur in the absence of an external magnetic 
field. Furthermore, weak SO interactions stabilize QSH 
states which preserve time-reversal symmetry. 

We considered the effects of long-range Coulomb 
and SO interactions. The omitted short-range part 
of Coulomb interaction can stabilize the TRS-broken 
spin-polarized states at / = 2. Taking the simplest 
Hubbard form of the short-range repulsion, Hjj = 
U^i=A B "iT"4i where the index i labels graphene lat- 
tice sites, we obtain a mean- field estimate of the en- 
ergy cost per particle (in the relevant pseudo LL) for 
QSH state. The result is 0.12[/[eV] meV, where we set 
B = lOT. Since the energy gain due to SO interactions 
has to exceed this value, for practical observation of the 
QSH state an artificial enhancement of the SO strength in 
graphene is likely to be necessary. This can be achieved, 
e.g., by graphene functionalization with adatoms [13, H^. 

The proposed states feature robust edge states (chi- 
ral in the case of AQH state, counter-propagating in 
the case of QSH state). We thus expect that the lat- 
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ter can be probed experimentally via the quantization 
of two-terminal conductance, much in the same wayit 
was done for the QSH state in HgTe quantum wells [13 ■ 
The interaction-induced gaps can also be probed in STM 
experiments. 
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K valley in the presence of the gauge field ([2]) can be 
rewritten in the following form: 



Hk — £0 



a 
at 



(20) 



where we have introduced raising and lowering operators 
of the magnetic oscillator [23|, which can be expressed 
via the complex coordinates z,z: 



A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 1097162 
(2009). 

E. -A. Kim, A. H. Castro Neto, Europhysics Letters 84, 
57007 (2008). 

S. V. lordanskii, A. E. Koshelev, JETP Letters 41, 574 
(1985). 

K. Sasaki, Y. Kawazoe, and R. Saito, Prog. Theor. Phys. 
113, 4637480 (2005). 

A. F. Morpurgo, and F. Guinea, Phys. Rev. Lett. 97, 
196804 (2006). 

5. V. Morozov et al., Phys. Rev. Lett. 97, 016801 (2006). 
V. M. Pereira, A. H. Castro Neto, Phys. Rev. Lett. 103, 
046801 (2009). 

Z. Ni et al., ACS Nano 3, 483 (2009). 

W. Bao et al.. Nature Nanotech. 4, 5627566 (2009). 

F. Guinea, M. L Katsnelson, A. K. Geim, Nature Phys. 

6, 30 (2010). 

F. Guinea, A. K. Geim, M. I. Katsnelson, K. S. 
Novoselov, Phys. Rev. B 81, 035408 (2010). 
N. Levy, S. A. Burke, K. L. Meaker, M. Panlasigui, A. 
Zettl, F. Guinea, A. H. Castro Neto, M. F. Crommie, 
Science 329, 544 (2010). 

N.-C. Yeh, M. L. Teague, S. Yeom, B. L. Standley D. A. 
Boyd, M. W. Bockrath, Surface Science 605, 1649 (2011). 
F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005) . 

K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 96, 
256602 (2006); M. O. Goerbig, R. Moessner, and B. 
Dougot, Phys. Rev. B 74, 161407 (2006). 
P. Ghaemi, J. Cayssol, D. N. Sheng, A. Vishwanath, 
[arXiv :1111.3640 (2011). 

B. A. Bernevig, T. L. Hughes, and S. C. Zhang, Science 
314, 1757 (2006). 

M. Konig et al, Science 318, 766 (2006). 
A. H. MacDonald, arXiv:cond-ma t/9410047] 

D. A. Abanin, D. A. Pesin, in preparation. 

Hongki Min, J.E. Hill, N.A. Sinitsyn, B.R. Sahu, Leonard 
Kleinman, A.H. MacDonald, Phys. Rev. B 74, 165310 

(2006) . 

A. H. Castro Neto, F. Guinea, Phys. Rev. Lett. 103, 
026804 (2009). 

C. Weeks, J. Hu, J. Alicea, M. Franz, R. Wu, 
larXiv:1104.3282. (2011). 

D. Ma, Z. Li, Z. Yang, Carbon 50, 297 (2012). 

M. Levin, A. Stern, Phys. Rev. Lett. 103, 196803 (2009). 
APPENDIX 



Landau level wave functions. Here we write down 
PLL wave functions, which were necessary for the anal- 
ysis of the interaction effects. The Hamiltonian of the 



a = -^{2£d + zl2i), ^ -^{-2id + z/2£). (21) 



V2' 



V2' 



Following Ref. [2y| , we introduce ladder operators &, , 
which commute with a, and therefore with the Hamil- 
tonian: 

b = -^i^C'd + z/2l), 6^ = ^("2^^ + (22) 

The eigenstates of the (non-relativistic) magnetic oscilla- 
tor are given by: 



n, m> 



V to! n! 



-|0), 



(23) 



where the vacuum state, annihilated by both a and 5, 



is given ipo.s)(z) 



e . Let us also introduce 



a notation (fn.m{z) = (zjn, m) for the \n,m) state wave 
function. 

The 71 ^ PLL wave functions of the Hamiltonian (|20|) 
are given by: 



^„,™(zj - V2 I sgn(n)^|„|,„(z) 



(24) 



The zeroth LL wave functions have a different form, with 
only bottom spinor component being non-zero: 



i'n..m{z) 







(25) 



The spectrum in the valley K' can be found in the 
same manner. One important difference is that the pseu- 
domagnetic field has an opposite sign, and therefore the 
operators a, b should be defined differently: 



7 ' 



^(2^9 + z/2£), flT = -y|(2''9 - z/2€), (26) 



b^ ^{2td + z/2l), 6^ = ^ 
^/2 V2 



{-2ld + z/2l). (27) 



The wave functions, defined using new operators a, b have 
the same form as those in the K valley p4|) (notice that 
this is different from the case of graphene in real magnetic 
field, where one would have to interchange the upper and 
lower spinor components). 



